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NUMERICAL  PREDICTION  OF  THERMAL  WAKES 


Introduction 

Wake  detection  has  long  been  considered  as  a  valuable  tool  for  the  location  and  identification  of 
ships,  i.e.  (1).  The  white  water  wake  was  first  used  for  this  purpose  (2)  by  aerial  observation  during 
the  second  world  war.  Later,  attempts  were  made  to  determine  if  the  thermal  disturbance  could  be 
detected  by  bow  mounted  thermistors  on  submarines  to  locate  surface  ships  (1).  More  recently  active 
and  passive  electromagnetic  remote  sensors  have  been  used  to  detect  the  various  elements  of  ship 
wakes  (3). 

Background 

Ships  moving  through  the  ocean  will  produce  disturbances  in  and  around  the  wake  which  can  be 
detected  using  a  variety  of  remote  sensors.  A  thorough  understanding  of  the  causes  and  the  evolution 
of  these  wake  disturbances  can,  in  principle,  provide  information  about  the  type  of  ship,  its  speed  and 
course,  and  probably  the  ambient  weather  conditions. 

The  first  major  ship-generated  disturbance  is  an  inviscid  Kelvin  wake  which  radiates  away  from 
the  ship.  This  wake  can  be  detected  either  visually  or  with  radar.  It  has  been  modeled  for  a  variety 
of  ships  and  ambient  sea  conditions  (4). 

The  second  major  type  of  disturbance  is  the  turbulent  wake  confined  to  a  narrow  region  behind 
the  ship.  Imbedded  in  the  turbulent  wake  are  vortices  created  by  the  bow  wave,  the  bilges  and  the 
propeller.  This  wake  has  been  the  subject  of  considerable  study  using  a  variety  of  techniques  includ¬ 
ing  inviscid  vortex  procedures,  linearized  Navier-Stokes  solvers,  and  more  recently,  a  parabolic 
Navier-Stokes  solver  with  a  k/epsilon  turbulence  model.  A  review  is  found  in  (3). 

A  third  disturbance  is  the  white  water  wake  caused  by  entrained  air  formed  by  the  bow  wave 
and  inside  the  turbulent  wake  by  the  bilges  and  the  propeller.  The  phenomena  of  air  entrainment  have 
not  been  modeled  although  the  evolution  of  a  bubble  wake  has  been  modeled  under  simplified  situa¬ 
tions  by  tracking  individual  bubbles  (6)  or  by  using  the  equations  governing  bubble  dynamics  (7). 

Another  type  of  wake  feature  is  the  thermal  ''scar'  formed  by  the  fluid  motion  inside  the  tur¬ 
bulent  wake  and,  to  a  lesser  extent,  changes  in  the  radiative  properties  of  the  free  surface  as  the  ship 
induced  current  interacts  with  the  ambient  wave  field  and  wind.  There  has  been  little  or  no  work  done 
to  model  this  wake  which  is  remotely  sensible  by  passive  radiometric  techniques  in  the  infrared  (IR) 
and  microwave  (MR)  bands. 

Observations  of  this  wake  show  two  distinct  parts  of  the  thermal  wake.  The  subsurface  wake  (1) 
is  fairly  predictable  and  remarkably  persistent,  often  lasting  60  minutes.  The  surface  scar  (3)  appears 
to  be  more  variable  and  is  less  persistent,  typically  disappearing  after  about  20  minutes. 
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A  more  complete  understanding  of  the  dynamics  of  turbulent  thermal  wakes  is  needed  as  an  ele¬ 
ment  of  a  detectability  assessment  for  surface  ship  wakes.  The  processes  by  which  thermal  wakes  are 
formed,  how  they  decay,  and  the  effect  of  ambient  wind,  temperature,  and  sun  are  not  well  known. 
This  report  contains  a  model  of  the  establishment  of  the  thermal  wake  and  evaluations  of  the  factors 
which  influence  its  development. 

Thermal  Structure  of  the  Open  Ocean 

The  variation  in  ocean  temperature  is  influenced  by  many  factors.  Even  without  the  presence  of 
external  disturbances,  such  as  ships,  the  ocean  temperature  is  rarely  uniform.  Typically,  the  variation 
of  temperature  with  depth  is  divided  into  several  regions.  The  upper-most  region  is  a  surface  micro 
,  layer  found  at  the  interface  between  the  ocean  and  the  atmosphere.  This  layer  is  typically  less  than  1 

mm  thick  and  is  characterized  by  the  dominance  of  molecular  interactions  over  turbulent  or  convec¬ 
tive  interactions.  At  this  laminar  layer,  heat  lost  to  the  atmosphere  through  long  wave  radiation  to  the 
sky,  convection  to  the  air,  and  latent  heat  loss  through  evaporation  is  balanced  by  conduction  to  the 
1  turbulent  mixed  layer  below. 

;  This  layer  is  also  associated  with  large  influences  of  surface-active,  monolayers  of  organic  films 

(8).  The  role  of  these  films  in  oceanographic  and  wake  processes  is  not  well  known  and  is  becoming 
a  subject  of  some  concern.  Film  effects  are  not  included  in  this  study. 

Below  the  surface  micro  layer  is  the  ‘mixed  layer’  which  is  dominated  by  turbulent  and  convec¬ 
tive  transport  properties  and  is  considered  to  be  uniform  in  temperature.  This  is  caused  by  the  rela¬ 
tively  vigorous  mixing  due  to  wind  and  wave  motion  at  the  ocean  surface,  called  thermocline  erosion, 

|  and  the  fact  that  the  large  majority  of  solar  radiation  is  absorbed  in  this  layer.  Below  this  layer  is  a 

thermocline  after  which  the  water  temperature  gradually  cools  until,  at  very  deep  locations,  the  water 
temperature  approaches  freezing.  The  mixed  layer  usually  extends  from  20  to  50  meters  below  the 
surface  and  will  be  the  deepest  region  of  interest  in  this  study. 

|  In  addition  to  variations  with  depth,  the  ocean  temperature  also  varies  in  the  horizontal  direction 

due  to  the  presence  of  oceanic  fronts  and  currents.  As  a  ship  passes  through  a  front,  the  characteris- 
|  tics  of  the  thermal  wake  can  change  dramatically  (3). 

If  the  mixed  layer  was  completely  uniform  in  temperature  there  would  be  no  observable  thermal 
I  wake.  The  bulk  mixing  due  to  the  passage  of  a  ship  would  cause  no  change  in  water  temperature. 

There  are  two  situations  where  an  appreciable  vertical  temperature  gradient  can  exist  in  the  upper 
mixed  layer. 

The  first  situation  occurs  when  wind  and  wave  mixing  is  weak.  Thermocline  erosion  can 
decrease  to  such  an  extent  that  solar  heating,  which  decays  exponentially  with  depth,  causes  an 
exponential  temperature  decay  with  water  depth.  This  temperature  variation  can  be  substantial  but  due 
to  the  extended  depth  of  the  profile  and  the  large  heat  capacity  of  water  it  usually  only  occurs  during 
the  afternoon  of  relatively  quiet,  cloudless  days.  During  the  day  solar  heating  effects  dominate  the 
infrared  signature  of  the  ambient  ocean.  The  mixed  layer  can  increase  in  temperature  by  several 
degrees  Celsius  and  reflected  and  scattered  solar  radiation  mask  emissions  by  the  surface  microlayer. 

The  second  situation  occurs  at  night.  Turbulent  and  convective  mixing  eliminate  temperature 
variations  in  the  mixed  layer.  Long  wave  thermal  radiation  produces  a  temperature  gradient  of  .1-.2 
|  C  in  the  surface  microlayer  where  thermal  conduction  dominates.  The  absence  of  solar  radiation  and 

temperature  variations  in  the  subsurface  mixed  layer  allows  this  temperature  profile  across  the  surface 
microlayer  to  become  easily  visible  in  the  infrared  spectrum. 
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Procedure 


As  mentioned  previously,  a  reasonable  model  of  the  turbulent  wake  behind  a  ship  currently 
exists  (9).  The  same  basic  approach  was  used  to  model  the  thermal  disturbance  behind  a  ship.  Instead 
of  the  finite  element  code  which  has  been  used  for  predicting  the  turbulent  wake,  a  finite  volume  code 
was  chosen  for  ease  of  code  development.  The  same  form  of  the  Navier-Stokes  equations  and  the 
Reynolds  stress  closure  was  employed  here.  The  parabolic  form  of  the  Navier-Stokes  equations,  the 
energy  equation,  the  turbulent  kinetic  energy  and  dissipation,  and  the  Reynolds  stress  closure  are 
shown  in  Table  1. 

This  form  of  the  parabolic  Navier-Stokes  equations  have  been  derived  in  (10).  They  are  formed 
using  the  usual  boundary  layer  ordering  of  terms  which  eliminates  the  axial  diffusion  of  momentum 
and  heat.  When  coupled  with  the  one-way  axial  flow  and  the  assumption  of  small  streamline  curva¬ 
ture,  the  elliptical  nature  of  the  equations  is  retained  only  in  the  transverse  directions. 

Integration  of  these  equations  begins  with  the  initial  velocity  and  temperature  field  behind  the 
ship.  All  equations  are  marched  downstream  using  time  (or  if  you  wish,  spatial)  steps  small  enough  to 
eliminate  the  need  for  iteration.  After  each  time  step,  the  elliptic  boundary  conditions  are  satisfied  for 
the  velocity  and  temperature  fields  only  in  the  transverse  direction.  The  continuity  equation  is  satisfied 
by  adjusting  the  transverse  pressure  distribution  using  a  Poisson  equation.  The  axial  pressure  distribu¬ 
tion  is  assumed  to  be  the  same  as  in  the  far  fluid,  in  this  case  zero  gradient. 


The  procedure  can  be  modified  when  the  effects  of  streamline  curvature  become  significant.  The 
three  dimensional  pressure  field  is  calculated  and  stored  during  each  sweep  through  the  calculation 
region  and  several  iterations  of  the  entire  calculation  are  required.  The  advantage  of  this  type  of  ‘par¬ 
tially  parabolic’  scheme  over  a  fully  elliptic  version  is  only  in  computer  storage  as  only  the  pressure 
is  stored  in  three  dimensions.  Execution  time  is  roughly  the  same  as  an  elliptical  procedure.  In  the 
case  of  ship  wakes,  it  may  be  possible  to  divide  the  computations  into  a  partially  parabolic  region 
near  the  ship  stern  and  a  fully  parabolic  region  extending  downstream  from  a  point  after  the  swirl 
velocities  become  sufficiently  small.  This  would  retain  most  of  the  computational  efficiency  of  a  para¬ 
bolic  scheme  while  improving  the  high  swirl  predictions  only  where  they  are  needed. 

The  energy  equation  is  slightly  simplified  when  compared  to  the  momentum  equations.  The  form 
of  the  energy  equation  suggested  in  (11)  uses  velocity-stagnation  enthalphy  correlations  to  model  the 
turbulent  transport  of  heat.  No  models  of  these  correlations  were  proposed.  Although  non-isotropic 
versions  of  these  correlations  can  be  found  elsewhere,  i.e.  (12),  it  is  unlikely  that  they  would  be  sig¬ 
nificantly  more  accurate  than  isotropic  versions  considering  the  different  modeling  assumptions  used. 
Therefore,  isotropic  fluctuating  velocity-temperature  correlations  were  included  in  the  present  calcula¬ 
tions. 


The  boundary  conditions  for  the  temperature  field  are  the  usual  cnes  (9)  except  at  the  free  sur¬ 
face.  There,  heat  is  added  or  removed  by  infrared  radiation  between  the  water  and  the  sky  (or 
clouds),  convective  heat  transfer  between  the  water  and  air,  and  latent  heat  transfer  proportional  to 
the  evaporation  rate.  All  three  effects  are  confined  to  the  surface  microlayer  of  the  water.  These 
effects  are  modeled  using  correlations  available  from  the  literature  (13).  Convective  sensible/latent 
heat  transfer  is  proportional  to  the  wind  speed  times  the  temperature/vapor  pressure  difference 
between  the  sea  surface  and  the  air  10  m  above.  Long  wave  thermal  radiation  is  constant. 


The  solar  flux  responsible  for  the  vertical  temperature  in  the  gradient  is  absorbed  over  a  much 
larger  region  of  the  ocean.  Depending  on  the  turbidity  of  the  water,  solar  heating  is  usually  restricted 
to  the  top  10  meters.  A  two  dimensional  simulation  of  the  ocean  surface  was  made  to  assess  the  time 
constant  characteristic  of  solar  heating  and  provide  initial  conditions  for  the  ship  wake.  The  transient 
temperature  profiles  shown  in  Fig.  1  indicates  that  several  hours  are  needed  to  influence  the  ocean 
temperature  significantly.  This  transient  lasts  much  longer  than  the  thermal  wake  is  expected  to  be 
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Table  1  —  Parabolic  Navier-Stokes  Equations,  Energy 
Equation  and  Reynolds  Stress 
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detectable  (1,3).  Therefore,  although  the  solar  influences  will  be  retained  as  source  terms,  their  effect 
on  the  calculated  thermal  wake  will  be  small. 

Additional  features  have  been  added  to  the  energy  equation  to  in  the  course  of  numerical  experi¬ 
ments  to  evaluate  the  effect  of  gross  anisotropy  of  the  thermal  diffusivity,  wind  stress,  buoyancy,  and 
the  nighttime  cold  layer.  These  modifications  will  be  described  later  in  the  report. 

Verification  of  the  code  was  made  using  ship  model  test  data  supplied  from  David  Taylor 
NSRDC.  Comparisons  were  made  of  the  axial  and  swirl  velocity  components  and  two  of  the  normal 
Reynolds  stresses.  The  agreement  was  reasonable.  No  temperature  measurements  were  available. 

A  benefit  arising  from  using  a  special  purpose  finite  volume  code  was  much  more  efficient 
storage  and  execution  time.  Initially,  the  Cray  XMP  supercomputer  was  to  be  used  for  this  study 
based  on  earlier  runs  made  with  the  finite  element  code  on  a  PDP  Vax  780.  The  control  volume  code 
ran  10  to  20  times  faster  than  the  finite  element  code.  Based  on  these  results,  the  code  could  be  easily 
used  on  a  minicomputer. 

Results  and  Discussion 

The  initial  velocity  conditions  used  for  this  series  of  thermal  wake  simulations  were  provided 
from  estimates  of  the  turbulent  wake  produced  by  the  Kaimalino,  a  small  waterplane  area/twin  hull 
(SWATH)  ship  operated  by  the  U.S.  Navy.  This  wake  is  slightly  simpler  than  conventional  ship 
wakes.  The  largest  part  of  the  hull  is  an  axisymmetric  body  which  is  completely  submerged.  This  is 
connected  to  the  superstructure  by  a  relatively  thin  strut.  A  horizontal  stabilizer  connecting  the  hulls 
and  a  foreward  canard  complete  the  underwater  structure. 

In  this  type  of  hull,  the  dominant  wake  structure  is  the  propeller  vortex  and  jet  (thrust)  superim¬ 
posed  on  the  drag  wake  of  the  axisymmetric  body.  This  wake  is  initially  completely  submerged  and 
only  interacts  with  the  free  surface  some  distance  behind  the  ship.  Small  necklace  vortices  will  be 
formed  at  the  joints  between  the  strut  and  control  surfaces  and  the  body.  A  weak  drag  wake  and  Kel¬ 
vin  wake  will  be  formed  by  the  strut  and  control  surfaces.  Due  to  limitations  of  the  empirical  scheme 
used  to  predict  the  initial  conditions,  the  necklace  vortices  and  Kelvin  wake  were  ignored.  The  tur¬ 
bulent  wakes  from  the  strut,  stabilizer,  and  canard  were  included. 

A  conventional  ship  hull  has  a  much  more  complex  shape  producing  strong  vortices  from  the 
bilges  which  interact  with  the  propeller  vortices  to  produce  a  multiple  vortex  pattern  behing  the  ship 
which  is  difficult  to  predict  using  empirical  methods.  In  addition,  the  Kelvin  wake  is  much  stronger 
and  transient  ship  motion  (roll  and  pitch)  can  also  produce  strong  wake  components.  Moreover,  all 
parts  of  a  conventional  ship  wake  are  formed  near  to  the  free  surface  which  means  that  not  only  are 
the  various  components  interacting  with  each  other  at  the  onset,  they  are  also  likely  to  ail  be  influ¬ 
enced  by  free  surface  conditions  and  processes. 

Thus,  using  the  SWATH  wake,  the  characteristics  of  the  thermal  wake  can  be  more  easily  traced 
to  particular  structures  of  the  turbulent  wake.  In  addition,  the  interaction  of  the  hull  with  the  free  sur¬ 
face  is  less  than  that  of  a  conventional  hull.  The  effects  of  the  interactions  between  the  free  surface 
and  turbulence  are  reduced  as  well  as  any  effect  of  the  Kelvin  wave  pattern. 

The  initial  velocity,  turbulent  kinetic  energy,  and  turbulent  dissipation  profiles  were  estimated  by 
superimposing  empirical  correlations  for  the  flow  past  simple  geometric  bodies  to  approximate  the 
underwater  portion  of  the  Kaimalino  hull.  The  initial  temperature  distribution  in  the  upper  ocean  was 
calculated  using  a  two  dimensional  control  volume  model.  This  model  assumes  constant  turbulence 
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parameters,  incident  solar  flux,  and  surface  heat  transfer.  Solar  radiation  absorption  as  a  function  of 
depth  was  calculated  using  low  turbidity  water  (14).  The  model  was  run  for  S  hours  of  solar  heating 
which  produced  approximately  a  1  degree  temperature  difference  between  the  water  surface  and  the 
layer  10  meters  deep.  This  profile  was  fitted  with  an  exponential  equation  and  then  used  as  initial 
conditions  in  the  model.  The  initial  isotherms  and  swirl  velocity  distribution  are  shown  in  Figs.  2a 
and  2b. 

Initial  Wake  Development 

Using  these  initial  conditions,  the  solar  source  term,  and  the  surface  heat  flux,  the  thermal  wake 
was  raieniaterl  starting  behind  the  ship  and  extending  3  km  downstream.  The  results  of  this  calcula¬ 
tion  are  shown  in  Figs.  3a  and  3b.  The  surface  thermal  wake  begins  some  distance  behind  the  ship 
when  the  propeller  vortex  rises  (or  expands)  enough  to  reach  the  surface.  Even  3  km  behind  the  ship, 
the  surface  wake  is  growing  in  intensity.  Although  not  shown  here,  this  behavior  continues  far  beyond 
the  3  km  station. 

The  effect  of  the  propeller  vortex  on  the  isotherms  is  shown  in  Fig.  3b.  The  swirl  velocity 
wraps  the  isotherms  into  a  spiral  and  the  effect  of  turbulent  diffusivity  is  enough  that  some  diffusion 
of  the  temperature  gradients  is  seen.  This  behavior  was  observed  qualitatively  in  (1). 

Infrared  observations  of  the  thermal  wakes  of  ships  (3)  show  the  wake  disappearing  within  3  to 
5  km.  In  situ  measurements  (1)  have  detected  subsurface  wakes  at  3  to  4  times  this  distance.  It 
appears  from  the  numerical  calculations  that  the  subsurface  wake  persists  because  the  primary 
mechanism  of  heat  transport  is  convective.  This  is  consistent  with  the  low  wind  speeds  necessary  for 
the  development  of  a  temperature  profile  which  causes  low  ambient  turbulence  intensities  in  the 
water.  In  the  absence  of  the  effects  of  a  strong  density  stratification,  the  propeller  vortex  which 
pumps  the  cold  water  to  the  surface  is  essentially  inviscid  and  decays  very  slowly.  The  surface  scar, 
however,  appears  to  have  another  mechanism  which  interferes  with  the  development  and  eventually 
causes  the  surface  scar  to  heal. 

Surface  Scar  Healing 

Several  mechanisms  which  might  account  for  the  eventual  healing  of  the  thermal  surface  scar 
were  identified.  A  numerical  scheme  was  developed  for  each  to  evaluate  its  potential.  The  results  of 
each  are  discussed  with  the  available  experimental  observations  in  the  following  paragraphs. 

From  the  initial  results,  it  appears  that  the  effective  thermal  diffusivity  has  a  small  role  in  the 
development  of  the  wake.  This  effective  diffusivity  is  modeled  using  the  turbulent  kinetic  energy  and 
dissipation  calculated  as  part  of  the  Reynolds  stress  model.  There  is  also  assumed  to  be  a  minimum 
background  diffusivity  in  the  ocean  (15).  Two  numerical  experiments  were  performed  to  estimate  the 
sensitivity  of  the  calculations  to  1)  the  background  turbulence  level  and  2)  the  anisotropy  which  may 
be  caused  by  shear  layers  i.e.,  wind  or  current.  The  interactions  between  the  free  surface  and  tur¬ 
bulent  anisotropy  were  not  modelled  here. 

The  first  experiment  involved  simply  increasing  the  minimum  value  of  the  thermal  diffusivity  by 
an  order  of  magnitude.  The  results,  shown  in  Figs.  4a  and  4b  show  the  isotherms  becoming  slightly 
diffused  over  a  larger  area.  Despite  the  larger  internal  heat  transfer,  the  scar  does  not  appear  to  heal 
more  quickly.  If  further  increased,  the  only  expected  result  would  be  a  smaller  temperature  gradient 
through  the  upper  ocean  layer. 
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The  second  mechanism  investigated  was  the  effect  of  anisotropy  in  the  thermal  diffusivity  caused 
by  horizontal  shear  layers.  This  approach  is  promising  because  preferential  heat  transfer  in  the  hor¬ 
izontal  direction  would  cause  the  surface  scar  to  heal  more  quickly.  Unfortunately,  it  is  also  expected 
to  cause  the  sub-surface  wake  to  heal.  This  is  not  in  agreement  with  the  limited  experimental  observa¬ 
tions.  The  results  of  a  two  order  of  magnitude  difference  in  the  horizontal  and  vertical  thermal  dif¬ 
fusivity  are  shown  in  Figs.  5a  and  5b.  Even  this  magnitude  anisotropy  causes  only  a  small  effect 
which  is  indeed  uniform  throughout  the  layer. 

Another  potential  mechanism  is  the  effect  of  perpendicular  wind  shear.  The  component  of  wind 
stress  perpendicular  to  the  ship  track  will  eventually  entrain  warm  water  from  areas  of  the  undisturbed 
ocean  upwind  of  the  wake  and  cover  the  surface  of  the  thermal  scar  leaving  the  sub-surface  less  dis¬ 
turbed.  While  this  explanation  appears  inviting,  it  would  only  be  a  factor  when  a  significant  com¬ 
ponent  of  the  wind  is  perpendicular  to  the  ship  track. 

This  effect  is  difficult  to  simulate  numerically  because  perpendicular  wind  stress  violates  the 
symmetry  condition  used  at  the  wake  centerline.  It  also  provides  an  inflow  boundary  condition  on  one 
edge  of  the  computational  domain  without  specifying  a  velocity  profile.  This  is  not  important  at  an 
outflow  boundary  because  information  is  being  convected  out  of  the  computational  area  unlike  the 
inflow  boundary  where  information  is  being  convected  into  the  domain. 

This  situation  was  modeled  by  removing  the  zero  normal  flow  boundary  condition  on  both  the 
vertical  transverse  plane  boundaries  and  substituting  a  zero  normal  velocity  gradient  condition  instead. 
The  alternative  was  to  extend  the  computational  domain  to  both  halves  of  the  wake  and  far  enough 
from  the  wake  so  that  a  logarithmic  velocity  profile  characteristic  of  drift  currents  could  be  imposed 
at  the  inflow  boundary.  A  computational  area  of  this  size  was  unreasonable.  The  effect  of  this  approx¬ 
imation  is  to  consider  only  the  turbulent  shear  forces  transmitted  vertically  from  the  wind  stress  at  the 
surface  and  ignore  the  momentum  convected  from  the  established  drift  current.  In  cases  where  the 
drift  current  is  either  very  deep  or  very  shallow  compared  to  the  wake,  this  will  be  a  good  assump¬ 
tion. 

Figures  6a  and  6b  show  the  predicted  effect  of  perpendicular  wind  stress  on  the  thermal  wake. 
Immediately  after  the  ship  the  wake  appears  and  begins  to  develop  as  with  the  standard  wake.  The 
wind  stress  begins  to  cover  the  wake  with  undisturbed  surface  water  at  about  2  km.  The  scar  com¬ 
pletely  heals  between  3  to  5  km  downstream.  The  combined  effect  of  the  swirl  and  establishing  drift 
current  on  the  isotherms  is  seen  in  Fig.  6b.  All  predictions  are  qualitatively  in  good  agreement  with 
both  sets  of  open  ocean  measurements  (1,  3).  One  set  of  informal  observations  made  at  the  David 
Taylor  Naval  Ship  Research  and  Development  Center  (NSRDC)  towing  tank  did  indicate  that  the 
wake  there  healed  eventually  without  the  assistance  of  wind.  There  are  other  differences  which  pre¬ 
clude  any  conclusions  based  on  this  observation.  The  ambient  turbulence  levels  were  not  representa¬ 
tive,  the  Reynolds  number  effects  are  very  different,  and  the  finite  width  of  the  towing  channel  itself 
could  be  a  factor.  Futher  experiments  there  could  prove  very  useful  in  validating  the  approaches  taken 
here. 
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Fig.  5a  —  Anisotropic  Thermal  Diffusivity  Isotherms 
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Fig.  6b  —  Surface  Temperature  with  Wind  Stress 


Buoyancy  Effects 


Based  on  observations  made  at  NSRDC,  a  mechanism  not  related  to  currents  or  wind  drift  can, 
in  some  instances,  cause  the  thermal  wake  to  disappear.  These  observations  were  made  in  water 
which  probably  has  low  levels  of  ambient  turbulence  and  measurements  taken  of  the  thermal  stratifi¬ 
cation  indicate  relatively  low  levels  of  natural  convection.  Because  all  forms  of  ambient  mixing 
appear  to  be  small,  a  flow  induced  mechanism  was  identified.  The  buoyancy  forces  acting  on  the  dis¬ 
placed  fluid  could  provide  such  a  mechanism.  In  some  cases,  density  variations  in  the  horizontal 
direction  can  cause  a  trailing  vortex  to  break  up  (16). 

The  proper  form  of  the  buoyancy  source  term  for  the  vertical  component  of  velocity  in  a  strati¬ 
fied  fluid  is  slightly  modified  from  that  found  in  a  fluid  of  uniform  density. 

5  =*  density  *g  *beta  *  ( T  -  T  (inf  ,z )) 

Thus  the  reference  temperature  is  not  a  constant  but  is  a  stationary  function  of  vertical  location  (17). 
This  form  of  the  buoyancy  force  has  been  used  only  in  semi-confined  flows  (18).  The  flow  pattern 
under  consideration  here  is  complicated  by  the  imbedded  vortex  which  can  become  unstable  under 
certain  conditions  (19). 

Simulations  using  the  initial  conditions  shown  in  Figs.  2a  and  2b  indicated  that  the  vortex  was 
destroyed  very  early  in  the  calculation.  There  was  little  of  the  deformation  of  the  isotherms  which 
were  predicted  earlier  and  which  have  been  observed.  The  kinetic  energy  of  the  transverse  velocities 
were  quickly  converted  to  potential  energy  by  the  small  displacement  of  the  stable  density  stratifica¬ 
tion. 

It  appears  that  this  initial  temperature  distribution  is  too  crude  to  account  for  the  effects  of  buoy¬ 
ancy.  A  first  order  estimate  of  the  effect  of  bucyancy  requires  a  well-mixed  vortex  core  temperature 
which  can  then  rotate  without  working  against  the  stable  density  stratification.  In  this  case  the  effect 
of  buoyancy  is  confined  to  the  fluid  entrained  by  the  vortex.  The  magnitude  of  the  effect  of  buoyancy 
will  depend  on  the  size  of  the  well  mixed  (isothermal)  vortex  core;  a  larger  isothermal  core  will  be 
less  affected.  Unfortunately,  there  is  essentially  no  information  on  the  initial  temperature  distributions 
in  the  early  wakes  of  wiips.  An  arbitrary  core  size  was  chosen  such  that  the  core  was  well  mixed 
where  the  swirl  velocity  was  .  1  m/s  or  greater. 

The  results  of  this  calculation  are  shown  in  Figs.  7a  and  7b.  The  surface  temperature  plot  shows 
a  reduced  signature  compared  to  the  standard  wake.  The  isotherms  at  5  km  show  that  the  effect  of  the 
vortex  has  been  reduced.  The  vortex  broke  down  early  in  the  calculation.  Energy  is  converted  from 
kinetic  rotational  energy  to  potential  buoyant  energy.  This  is  in  addition  to  the  effect  the  entire  vortex 
structure  rising  or  descending  due  to  buoyant  forces. 

The  healing  mechanism  is  seen  starting  at  the  end  of  the  computation.  Warm  water  is  being 
drawn  over  the  scar.  Even  after  the  primary  scar  is  covered,  a  less  well  defined,  wider  scar  will 
remain.  Predictions  based  on  this  calculation  qualitatively  fit  the  available  ocean  measurements  (3). 
Quantitative  comparisons  are  not  possible  because  of  the  differences  between  the  ships  used. 

Night  Wake 

A  ship  moving  through  the  night  ocean  produces  a  completely  different  wake  than  during  the 
day.  Without  the  solar  heating,  the  largest  source  of  heat  transfer  is  removed.  Ambient  turbulent  mix¬ 
ing  quickly  eliminates  the  temperature  profile  in  the  upper  ocean.  Heat  transfer  to  the  air  occurs  only 
in  the  top  few  microns  of  water  primarily  through  re-radiation  to  the  night  sky,  latent  heat  lost  by 
evaporation,  and  convective  heat  transfer  with  the  air.  The  relative  importance  of  each  of  these  terms 
depends  on  the  atmospheric  and  ocean  conditions.  Usually  a  cooled  layer  forms  on  the  top  of  the 
ocean.  Such  a  cool  layer  is  unstable  except  at  the  interface  between  the  water  and  air  where  turbulent 


fluctuations  are  damped  and  a  viscous/conductive  region  forms.  As  the  cool  layer  tries  to  grow 
beyond  this  layer  Rayleigh  type  instabilities  reduce  the  layer  to  the  viscous  thickness. 

This  layer  has  been  observed  in  the  laboratory  and  open  ocean  (13,20).  It  remains  intact  through 
a  variety  of  weather  conditions  and  apparently  breaks  down  only  at  the  onset  of  wave  breaking  or  a 
similar  vigorous  surface  disturbance  such  as  a  ship  passage. 

This  wake  was  modeled  here  by  reducing  the  surface  control  volumes  to  the  observed  thickness 
of  the  viscous  layer.  A  vortex  will  pump  water  to  the  surface  causing  the  viscous  layer  to  thin  until  it 
reaches  the  temperature  of  the  subsurface  water.  The  effect  will  be  a  warm  wake  on  a  cold  surface. 
The  thermal  scar  would  heal  when  the  vortex  weakened  enough  to  permit  heat  transfer  to  the  air  to 
again  dominate  the  thinning  effect.  Because  the  water  is  at  uniform  temperature  except  at  the  top 
node,  the  energy  equation  is  solved  only  for  the  surface  nodes. 

Figure  8  shows  the  predicted  form  of  the  thermal  wake.  It  appears  to  be  reasonable  when  com¬ 
pared  to  the  rare  (and  undocumented)  observations  (21).  Laboratory  studies  (13)  indicate  that  this 
mechanism  (dilution)  is  not  the  cause  of  the  extended  wake.  The  only  known  mechanism  capable  of 
causing  the  disruption  of  the  cold  layer  is  the  surface  breaking  (breaking  waves,  the  passage  of  a 
ship).  When  these  mechanisms  cease,  the  wake  re-establishes  within  20  seconds.  Obviously,  the  thin¬ 
ning  of  the  layer  alone  is  insufficient  to  explain  the  observations  of  the  night  wake  phenomenon. 

Two  possible  mechanisms  which  are  currently  being  investigated  involve  the  special 
biological/chemical  characteristics  of  the  surface  microlayer  (e.g.  see  ref.  20)  and  the  disruptive 
behavior  of  the  bubble  wake  on  the  sea  surface.  The  biological  and  chemical  surfactants  can  accumu¬ 
late  on  the  surface  of  the  sea  due  to  the  ship  induced  vortex  and  the  motion  of  the  bubble  wake.  Sur¬ 
factants  accumulate  preferentially  on  the  air/water  interface  which  makes  bubbles  particularly  effective 
of  concentrating  surfactants.  Once  on  the  sea  surface,  the  surfactants  would  disrupt  the  surface 
micrlayer  until  natural  diffusion  processes  could  disperse  the  slick.  r  bbles  can  also  distrupt  the  sur¬ 
face  microlayer  mechanically  by  rising  to  the  surface  and  bursting. 

Summary  and  Conclusions 

The  numerical  model  developed  here  represents  a  fully  parabolic  description  of  the  turbulent 
wake  with  anisotropic  Reynolds  stress  closure  and  isotropic  turbulent  energy  closure.  Predictions 
indicate  that  the  initial  thermal  wake  development  is  in  reasonable  agreement  with  the  limited  observa¬ 
tions.  The  mechanisms  which  cause  the  eventual  healing  of  the  surface  scar  while  the  subsurface  ther¬ 
mal  disturbance  persists  were  investigated  using  numerical  experimentation. 

The  most  promising  physical  mechanisms  for  the  observed  wake  behavior  are  the  effect  of  per¬ 
pendicular  wind  stress  and  buoyancy.  Wind  stress  appears  to  be  a  good  candidate  for  thermal  wake 
healing  whenever  a  significant  component  of  the  wind  is  normal  to  the  wake.  The  effect  of  buoyancy 
was  difficult  to  predict  quantitatively  because  of  vortex  breakdown  which  appears  to  occur  early  in  the 
wake  development.  Again,  this  result  is  preliminary  because  good  data  for  comparison  are  not 
currently  available. 

The  night  wake  was  predicted  assuming  that  surface  spreading  caused  the  breakdown  in  the 
viscous  cold  layer.  Although  this  assumption  is  debatable,  the  predictions  appeared  reasonable  com¬ 
pared  to  limited  and  qualitative  descriptions  of  the  phenomenon. 

For  both  day  and  night  wakes,  good  quantitative  experimental  data  of  the  surface  and  subsurface 
thermal  wake  is  necessary  to  verify  these  predictions  or  to  provide  a  basis  for  further  model  develop¬ 
ment.  In  addition,  a  more  basic  understanding  of  vortex  propagation  in  a  stratified  fluid  is  required  to 
establish  a  relationship  between  buoyancy  and  vortex  breakdown.  Finally,  the  basic  dynamics  of  cold 
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surface  layer  breakdown  and  re-formation  is  also  needed  for  guidance  in  developing  more  complete 
physical  models  of  the  night  wakes  of  ships.  In  view  of  the  scaling  difficulties  and  other  limitations 
associated  with  laboratory  experiments  it  may  be  necessary  to  pursue  experiments  at  or  near  full  scale 
under  carefully  selected  and  documented  conditions. 
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